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ABSTRACT 


A 

* 


A general Fokker-Planck equation is deduced, which describes the distribu- 
tion of geomagnetically trapped electrons as a function of longitude, time, energy 
and mirror point field intensity. A special variable for the longitudinal position 
of a particle is introduced. The coefficients representing longitudinal drift, ioni- 
zation loss and multiple Coulomb scattering in the Fokker-Planck equation are 
derived. The physical interpretation of this equation is analyzed for several 
special cases. It is found that the usual procedure of averaging over longitude 
in order to obtain a longitude — independent description is not valid for electrons 
mirroring at low altitudes in the South American Anomaly. The correct procedure 
is discussed. It is further concluded, that the region East of the Anomaly, ini- 
tially depleted by precipitation, is replenished by electrons whose mirror points 
were situated in a narrow "window" in B, before passing through the Anomaly. 
A very limited extension of the atmosphere in the Anomaly should control this 
mechanism of replenishment by Coulomb scattering. A small diurnal effect for 
the low altitude electron flux is predicted, determined by the diurnal variation of 
the atmosphere in the Anomaly. 
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LONGITUDE DEPENDENCE OF GEOMAGNETICAL LY 


TRAPPED ELECTRONS 
PART I 


I. INTRODUCTION 

Time has come to look with more detail into the latitude dependence of geo- 
magnetically trapped particle fluxes from the theoretical point of view. Until 
recently, experimental information on this subject was very scarce. Only after 
the discovery of a region over the East coast of South America and the adjacent 
Atlantic Ocean, with enhanced fluxes of trapped radiation at a given altitude 
(Yoshida et al. 1 , Vernov et al. 2 ), attention got more and more focussed on the 
experimental analysis of particle precipitation in the Anomaly, and the mechan- 
ism of subsequent replenishment of the depleted shell regions. In particular, 
Imhof and Smith 3 , Paulikas and Freden 4 and Mihalov et al . 5 have made a care- 
ful study of artifically injected electron fluxes, measured by various satellites. 
On the other hand, Freden and Paulikas 6 and Heckman and Nakano 7 have ana- 
lyzed proton fluxes at low altitudes in the Anomaly. Evidence for Bremsstrahlung 
X-rays from electrons precipitating into the Anomaly was found by Ghielmetti, 
et. al. 8 . A thorough study of the longitude dependence of electron fluxes was 

9 

published by Williams and Kohl . It is therefore desirable to set up a theoreti- 
cal description of the longitudinal behavior of trapped particles, and to test by 
comparison with experimental data, the various assumptions made about inter- 
action processes governing particle diffusion. 

Theoretical description of electron trapping, diffusion and precipitation was 
so far done only for configurations averaged over all longitudes 10 - n >12,13,14. 
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In these papers, a time-dependent Fokker- Planck equation was set up for the 
electron distribution function and used to follow the evolution in time of a given, 
initial electron flux. 

In order to study the longitude dependence of the electron distribution on a 
given magnetic shell, it is necessary to derive a more general diffusion equation 
which contains an additional variable related to longitude. In this equation the 
coefficients, which depend upon the atmosphere, will be a function of longitude 
and local time of day in the Anomaly. One important by-product of the equation 
derivation will be the identification of an appropriate variable by which to de- 
scribe the longitudinal dependence. We will call this variable X. It will be used, 
together with the particle energy E, the scalar magnetic field B at the mirror 
point, and the well known Mcllwain shell parameter L, ,s to describe the space 
in which particle densities change with time, t. 

The main purpose of part I is to set up such a general equation, to discuss 
its physical meaning, to compare it with the previously used longitude-averaged 
equation, and to draw some general, qualitative conclusions. The second part 
will deal with atmosphere-field configurations to be used in this longitude de- 
pendent description, and results of a numerical integration of the general equa- 
tion will be presented. 

Before setting up our equation, let us picture the problem in general terms. 
Consider the familiar B-L space, 15 in which trapped radiation fluxes are usually 
described (Fig. 1). It can be shown that electrons whose mirror points are be- 
low 100 km cannot remain trapped for more than a few bounces. Thus, it is 
qualitatively useful to consider the 100 km level as the location of a sink. We 
shall return to this point in more detail later. At a given longitude X the 100 km 


1-2 




Figure 1— The “windshield wiper” effect in B-L space. 

level can be displayed as a locus B c (L,X) in B-L space. Essentially no electrons 
will be found with mirror points above this curve, that is, for mirror point fields 
larger than B c (L,X\ We have plotted two extreme loci of the 100 km level, cor- 
responding to a longitude (that is, a value of X) right in the "center" of the South 
American Anomaly, and a longitude over the Pacific Ocean, respectively. From 
Fig. 1, one can see that for particles of a given energy, the 100 km curve, 

B c (1, X), "oscillates" with longitudinal drift frequency between these two extreme 
positions. As B c (L,X) lowers from its maximum, or "Pacific," position, the 
region in B-L space between the two extreme positions, called the "shadow re- 
gion," is wiped clean of particles. As B c (L,X) rises from its minimum, or 
"Anomaly," position, the opportunity exists for atmospheric scattering, or any 
other nonadiabatic process, to repopulate this shadow region. This wiping clean 
and repopulation has been likened to the action of a windshield wiper. 
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Satellite observations do, in fact, find electrons in the shadow region 3,9. 
This observation means that particles have enough non-adiabatic interaction 
during one drift around the earth to diffuse into the wiped-out or "shadow" 
region. Our main goal is to study this process of fast replenishment 
and to determine to what extent it can be explained by the process of in- 
teraction of electrons with the atmosphere through Coulomb scattering and 
energy loss. 

II. REPRESENTATION OF THE ELECTRON DISTRIBUTION 
AS A FUNCTION OF LONGITUDE 

We shall describe the trapped electrons by the distribution function used by 
Welch, Kaufmann and Hess u . We shall denote by 

SN = U(I, B, X, E, t) 8<D8B8E (1) 


the number of electrons contained at the time t in a tube of field lines of mag- 
netic flux 8$ = B Q SA , situated at a longitude characterized by the (not yet de- 
fined) parameter X , with mirror points between B and B+SB, and with energies 
between E and E +8E (See Fig. 2). In this paper, the mirror point field will al- 
ways be denoted by B without suffix. The shell on which these electrons drift 
is characterized by the two adiabatic invariants /j- and 3 16 , which for static 
fields and in absence of an electric field reduce to the field-geometric invariants 


B 


(mi rror ) 


const. 


(2a) 


9 



(2b) 
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respectively. For the time being, we shall make no limiting assumption regard- 
ing azimuthal symmetry; this means that particles initially on one and the same 
field line may not at all stay on a common shell during their drift in longitude 17 ' 18 • 



The distribution U defined in (1) is in general not directly measurable. Rela- 
tions with other particle distribution quantities are given in u . In particular, U 
is related to the directional, differential flux at the equator j Q = j (I, a 0 , X, E, t) 
by the expression 


ii 2?t • 

U = — Jo T : 


B 2 


0 b 


( 3 ) 


a o = arc sin Vb o /B is the equatorial pitch angle of a particle mirroring at B; r b 
is the bounce half-period (time to travel from one mirror point to its conjugate). 
In an azimuthally asymmetric field, the pitch angle a Q of a given particle is in 
general a function of longitude a- 0 - a Q ( B, X) (in these fields, the field intensity 
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B 0 at the equatorial ring of a shell varies with longitude). This is why we pre- 
fer to use the distribution function U, which is represented in terms of the longi- 
tudinally invariant variable B. 

Let us now consider electrons of a given energy, mirroring at the same B- 
value, trapped between two neighboring shells labeled with the values of the 
second invariants, I and I +5 I, respectively. I+SI is the value of the second 
invariant, taken along a line of force belonging to the outer shell, between mirror 
points with the same field intensity B. We shall call V Q i = n S i/Sy the gradient 
of I at the equatorial point of a magnetic shell (Figs. 3, 4). Northrop and Teller 16 
have shown that the equatorial drift velocity, averaged over one bounce, is 
given by 


u o r 


mv _ _ p V 0 I 

— — V 0 IxB= b 

eB o r b eB o T b 


where p is the particle's momentum, e the charge (<0 for electrons) and r b 
is the half-period of bouncing. B Q is the equatorial B -value. For the different 
vectors, see Fig. 4. Introducing S b , rectified path of the particle between mir- 
ror points, we can write 

m o c2 V Q I 

u 0 =/3 2 7 b (4) 

eB 0 S b 

fi and y are the usual relativistic factors. Notice the following useful general 
expression for the path S b , which can be obtained straight-forward by taking the 
derivative of (2b) with respect to the mirror point field R, along a given field 
line: 

^ = I + 2B |* (5) 


1-6 



We denote with Sx = u Q St the distance between equatorial points of two 
neighboring field lines of a given shell; Sx is then the element of arc of the 
shell's equatorial B Q -ring; x = fdx is the total arc length along the equator 
from a given initial point and can be used as a label to locate field lines on a 
given shell. Evidently 


x 



dx 

— dy 

OCp 


(5) 


where Bx/Bcp is a field-geometric factor (3x/9<p = 27rRL/360 for a centered 
dipole) , and 9 the geographic longitude of the equatorial point. Notice that for 
the angular equatorial drift velocity, 


<Po 


dx/dy 


u , 


( 6 ) 


Let us take a tube of lines of force of equatorial cross section SxSy and mag- 
netic flux B 0 SxSy (Fig. 3). By the definition of U, (1), the number of elec- 
trons in this flux tube, which mirror between B and B+SB and with kinetic 
energy between E and E+SE , is given by 

SN = UB 0 S x 8y SBSE (7) 

We now want to follow the history of these particles, as a function of time. 
Before setting up the Boltzmann equation governing the distribution U , we have 
to introduce a convenient variable for the longitude, i.e., a convenient label for 
the field line of a given magnetic shell, around which a particle is instan- 
taneously spiralling. We cannot take y or x (5) as suitable variables, because 
U is a distribution in flux , not in <p or x. And in the flux expression B 0 8 x 8 y 
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intervening in (7), not only Sx, but also Sy is a function of longitude, for the 
general geomagnetic field. 

In order to find the correct longitudinal variable, we have to transform the 
flux into an expression SO = SXSY, in which SY - still related to the radial in- 
crement Sy — is now longitude-independent, so that it may be ignored when fol- 
lowing the particles during their longitudinal drift. In that case, X would be the 
correct longitudinal variable, its differential SX containing complete information 
about the longitpde dependence of the flux SO of a tube filled with particles, as 
they drift around the earth. It is easy to realize that, in our notation, X and Y 
must be proportional to a pair of " a,/3 " coordinates 19 . 

In order to find X and Y, let us write the obvious relation (Figs. 3 and 4) 


Introducing this into the expression for the flux element, we have: 

SO = B 0 Sx Sy = B 0 Sx^; 

v o i 

A possible pair of coordinates is therefore 

X= / B »v < 8a > 

Y = I (8b) 

The integral (8a) is to be performed along the equatorial ring of a given I , B - 
shell, starting at a fixed reference longitude. We must emphasize again that B Q 
is in general longitude dependent, and therefore cannot be taken out of the integral 
(8a). Notice further that our system X, Y is independent of the particle's energy, 
depending only on the field geometry. 
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4 

Using (8a) and (8b) , the flux element becomes 

§4> = B q SxSy = SX8Y 


where SY (=81) is now independent of longitude and can be left out in the entire 
following discussion. With this new variable X , the longitudinal drift velocity 
becomes, according to (8a) and (4): 


i2L = u B o = _p ^R 2 y J, 

? dt U °dx U °V^ eT b e P 7 s b 


( 9 ) 


Remember that S b is in general longitude dependent. Notice the following rela- 
tion, taking into account (5) and (6): 


| T b dX = ^ L= T 2 ( 1Q ) 

0 9 0 

m. THE GENERAL FOKKER-PLANCK EQUATION FOR TRAPPED ELEC- 
TRON DIFFUSION 

We are now in condition to set up the general equation governing the distri- 
bution of trapped electrons. These electrons will undergo displacements in B , 

E and longitude, caused by three types of mutually independent interactions: 

(la) A change in I, coupled to a change in B due to the stochastic process 
of pitch angle scattering (multiple Coulomb scattering in the atmos- 
phere, interactions with traveling electromagnetic waves, etc.); 

(lb) A change in I , B and E due to the stochastic process of radial diffusion 
across shells (violation of the third adiabatic invariant during geomag- 
netic perturbations, etc.); 

(2) A change in kinetic energy E, due to ionization loss in the atmosphere, 
considered here as a non - stochastic process (i.e., neglecting straggling); 
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(3) A change in longitude due to interaction with the static magnetic field 
(longitudinal drift), again a non-stochastic process. 

All these processes are physically in dependent of each other, although all inter- 
vening parameters are in general functions of the three variables B, X and E. 

Let us relate the distribution of electrons as it appears in (20), with the 
distribution of electrons at a slightly earlier time t -At , at a different longi- 
tudinal position. For the time being, we shall forget the non- stochastic charac- 
ter of E and X, and treat all variables as if they were of the same, stochastic 
nature. We write: 

U(I, B, X, E, t) = 

JJjJ U ( 1 - T 7 , B-/3, X-£, E-e, t-At)Il (I- 17 , B-/3, X-£, E-e, ( 11 ) 

t), /3, g, e, At)dr?d/3d£de + Q (I , B, X, E, t)At 

In this relation, the distribution of electrons which at the time t are at a posi- 
tion X , is linked to the distribution of those electrons which at an earlier time 
t - At were at X - £ , and which happened to diffuse, scatter, slow down and 
drift the right amount in the interval At, in order to become part of the 
population described by the left hand of ( 11 ). 

Fldryd/Sd^de is the a priori probability that these electrons have undergone 
just the right changes in I, mirror point field, longitudinal position and energy, 
in the time interval At . 

The source term QAt represents the contributions of electrons added by 
injection to the original bunch of particles, during At. 
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Expanding- all intervening functions in Taylor series in ??, fi 9 £ , e and At, 
we obtain: 


s? + a <u ^>> = n <0 ^ )t i^ <u ^> : 


+ cross terms and higher order terms 


The brackets stand for the average per unit time of the enclosed variables: 

= 5iJIfJ' nd,d/3dfd £ 

We can now re-instate the non-stochastic character to X , by taking n as a 
delta function in longitude: 

n ( I , B, X, E, 77, /3, <f, e , At) = U Q (I, B, X, E, tj, J3, e, At) 8 (£ ~t At) 

€ is given by (9). n Q is the probability for a change in I, B and E due to radial 
and pitch angle diffusion, in the interval At. It is conveniently expressed by 

n (I, B, X, E, 7 ), /3, e, At) = P ^ 1, B ’ X ’ E ’ 1 1 ) A t (12) 

0 T b ( I » B, X, E) 

where P is the scattering probability per one half-bounce . In all this we have 
implicitly assumed that scattering effects are extremely small during one 
half-bounce of the particle. In other words, we suppose that P« 1 for -q, / 3, e ,1 0 
and approximately P = 1 for 17 , J3, e , = 0. Equation (11) finally becomes: 

"4 <U ^ >+ 23B^ <U ^^ ) 'm <U ^' >) + 0+ <13> 

+ cross terms and higher order terms 
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This is the most general equation of diffusion for trapped particles, applicable 
for an arbitrary trapping field geometry and for arbitrary radial and pitch angle 
diffusion processes. All coefficients ( ) (average changes per unit time), can be 
expressed as average changes per half-bounce, in the form: 

(* 1 • ) = 7- < > (I 4 ) 

b 

where 


{ } = 



• • Pd7jd/3de 


( 15 ) 


These coefficients are not independent of each other. For instance, if we con- 
sider the action of a radial (cross-shell) diffusion process, any independent 
change of, say, I , will be accompanied by dependent changes of mirror point 
field intensity and energy, determined by the condition of conservation of the 
first and second invariants p. and 3 . On the other hand, for a pitch angle 
scattering mechanism, any independent change SB of B will be accompanied by 
a dependent change Si of I, determined by the condition for the particle to re- 
main on the same field line during the scattering process (within one gyroradius, 
of course) : 


SI = — i— SB 

Bl/BB 

The derivative B I/BB is to be taken along the field line at which the scattering 
occurs. This leads to the relations 


<**>-( urn-- 



for any pitch angle scattering process. This allows a considerable simplifica- 
tion for nearly azimuthally symmetric fields. Let us change variables in the 
form L = L(I, B), B' = B, X' -X. We can rewrite the Fokker Planck equation in 
these variables, introducing new coefficients, in the form: 


f * aT <u|) = - k <“<*> > + 5 £ < u <^> - si > + 


3B 


etc. 


The relation between (k) and (j3) for pitch angle scattering processes will 
now be of the form: 


(*) = (^ 2 ) =(|) 2 <^ 2 >' etc - 
Again, the derivatives are taken along the field line. For a given field line, it 
is always possible to find a relation L = L(l , B) such that dh/dB - 0 for all 
points of this line. This condition, ingeneral, will not be fulfilled for other lines 
of the same shell. However, for nearly azimuthally symmetric fields, it is 
possible to find a unique function L such that 3L/3B = 0 for all field lines of a 
given shell. This condition implies degeneracy of all shells passing through 
one field line 17 (particles mirroring on one field line, populate the same surface 
during their drift motion, irrespective of their pitch angles or mirror points). 
For the (internal) geomagnetic field, this function is Mcllwain's L parameter. 

Using the L -parameter representation instead of I, the corresponding co- 
efficients (k) , (k^ , etc. vanish in the Fokker-Planck equation, for pitch angle 
scattering in nearly azimuthally symmetric fields. As there is good evidence 
that the geomagnetic field does have such properties within about 4-5 earth 
radii, we shall use this representation in all what follows, and therefore drop 
the parameter L from our expressions. Before writing our Fokker-Planck 
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J 


equation in its final form, we shall re-instate the non-stochastic character to 
the energy, assuming that the only process for energy change is that of ioniza- 
tion loss, without straggling. This means that the probability P in (12) can be 
expressed as a delta in energy too: 

P(B, X, E, A e) = P'(B, X, E, /3) 8 (E - e At) 


where e is the energy loss per unit time (e < 0 ). 

The final form of our Fokker-Planck equation is therefore, up to the second 

order term: 

^ t li <u#) = -i (U ^ ) + 2^ (U ^ ) -| <ui) + Q <16) 

We can re-write the longitude-convection term as a function of ordinary 
geographic longitude, 9 , taking into account (9) and ( 6 ): 




Notice that this is not equal to V ^9 ( 9 0 U), as one would have erroneously 
obtained, had one written the Fokker-Planck equation in the variable 9 . 


(17) 


IV. DISCUSSION OF SPECIAL CASES 

Let us now discuss equation (14) from the point of view of its physical 
meaning. First of all, in absence of any diffusion processes and sources, we 
have: 


By 

Bt 


+ 


^<ui) = o = flU<p oTb ^(u/ Tb> 


which can be written 


_3_ /U_ 
Bt \r h 


0 B 9 


- 2 - 1=0 


(18) 
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Remembering that U/r b is proportional to the directional flux j 0 (3), the above 
expression represents Liouville's Theorem: a flux j 0 , injected at any longitude, 
proceeds drifting along the shell, maintaining its value unchanged, provided one 
looks into the correct direction, i.e., along the actual path of the given class of 
particles. For the earth's field, equation (18) represents the "flight-time spec- 
trometer" effect for a group of particles artificially injected at high altitudes. 

If we now integrate equation (14) over one complete longitudinal cycle of X 
for a fixed time t, the second term on the left vanishes. We obtain: 


Tt = f u ^> dx + i^u<S’> d x--£ju*dX + Q 

dividing by <j) dX and calling 

U . „ - *—2 


dX 


dX 


f 

feu 

e AV =- L jr- 

r 


dX 

dX 


dX 
U dX 




A V' 


| (/3 2 ) U 

j)U dX 


dX 


We obtain the longitude-independent equation 

B 


w *' = -£-<« 


at 3E '~ AV c Av) "bb (U av(^)av) + 2 Bi~2 ( u av (^ 2 )av> + Q 


( 19 ) 


This equation is formally equivalent to the time -dependent equation used by 
several authors 10 ' 11,12,14 . Notice, however, the following remarks: 
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(1) The distribution U AV is averaged over the new variable X. 

(2) The "coefficients" ^ AV ,^) AV and ^ are not simple werages 

over the new coordinate, but they are weighted with the distribution U 
itself. In other words, they are functionals of the unknown distribution. 

All this leads us to the conclusion that a longitude-average Fokker-Plarck 
treatment of the problem of trapped particle diffusion is not valid, unless we re- 
strict our description to only those particles which mirror high enough at all 
longitudes so that thei r distribution function can be expected a priori as very 
little longitude dependent. Equation (19) is certainly meaningless as a differ- 
ential equation for that portion of B -L space, which descends below about 200 km 
in the South American Anomaly. This is precisely the domain where the longi- 
tude-independent treatment has so far failed to give numerical results compatible 
with experimental measurements on particle loss rates. 

In order to "legalize" the longitude-independent description for such cases 
in which we know a priori that U will not depend strongly on X , we have to 
evaluate the coefficients of (19) for a distribution U which has the smallest pos- 
sible longitude dependence. This is given by the solution of (18;: U/r b ~ j 0 - 
const. Taking into account (9), we obtain for the coefficients of the longitude- 
independent equation: 



1-17 


These expressions are formally equivalent to the longitude-average coefficients 
used by Hassitt 20 . They represent a longitude average, weighted with the inverse 
of the local drift velocity. One important difference should, however, be noted: 
due to the longitude-dependent asymmetry of the real geomagnetic field, integrals 
along the field line contained in e ^ /3) and cannot be permuted with the 

integration over longitude in (20). A correct evaluation of the longitude-inde- 
pendent coefficients therefore requires first, separate evaluation of each co- 
efficient for a certain number of shell field lines, and only then, weighted aver- 
aging over longitude. 

We must finally insist again that the use of a longitude-independent Fokker- 
Planck equation for the description of trapped electron diffusion, is not valid at 
all for those high B values, for which the atmospheric interactions in the region 
of the Anomaly will cause a notable longitude dependence of U, or j Q . 

We shall now discuss equation (16) for a steady state 3U/Bt = 0, as it is the 
case for natural inner belt electrons during geomagnetically quiet epochs. We 
then have 


_3_ 

BX 


(U£) 


,9 


3<p 



<U</3’>>-=|< u <VQ (21) 


This will be the basic equation for all what follows: The solution of (21) is sub- 
ject to a natural boundary condition at the equatorial B-value: 


3U 

3B 


= 0 for B = B 0 


( 22 ) 


At the ends of the field lines in the high atmospheric density region, there is no 
physical boundary condition at all: as particles enter higher densities, they get 
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slowed down by ionization loss and die away in energy space. However, for 
practical reasons, one can impose an artificial boundary condition 

U = 0 for B > B c (X) (23) 

which represents a sink of electrons at a level where the field intensity is B c 
(usually taken as the B value at 100 km altitude). 

The solution of equation (21), with boundary conditions (22) and (23), de- 
pends on the source term 0. Unfortunately, very little is known about this term, 
for inner belt electrons. For sure, its contribution is very small during time 
intervals of the order of one longitudinal drift period. But this is just the 
typical time scale of importance when one wants to study the longitude de- 
pendence of the electron flux. Of course, one has to start with a given distri- 
bution which must be close to an eigenmode of the general equation (16). For 
the time being we shall neglect the source term in front of all other terms of 
equation (21). ). 

A useful concept for a qualitative analysis of equation (21) is that of "mirror 
point flow" along a field line. Let us consider again a population of trapped 
electrons, mirroring between B and B + SB, enclosed between two surfaces 
labeled I and I + SI (Fig. 3). According to (1), the number of electrons, drifting 
through a given longitude (labeled X) during a time interval St, will be given by 

U — SYSBSESt =UcfSBSESY 
dt 

The total number of electrons mirroring between B and the equator, drifting 
through a given longitude per unit time will then be 
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f B . 

N= U£dB 

K 

(We have divided by the constant increments SE8Y). We shall define as "char- 
acteristic mirror point trajectory," the function B m =B m (X). such that the above 
integral remains constant throughout the drift motion (Fig. 5): 

r B « (x) • 

N= I U^f dB = const. (24) 



B m (X) of course does not represent the actual path of a particle's mirror point. 
However, if one neglects energy loss and source contributions, it becomes the 
average path of the mirror point of a particle in B-X space. This even is true 
if one does take into account energy loss; only that in this case particle individu- 
ality would be lost: a given particle lost out of the energy interval, might be re- 
placed by another one flowing in from a higher energy. Making use of (24), and 
(10), we now introduce the quantity 
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£ (UO dB 


(25) 


dB dB dX 

on m 

dcp dX dtp 


- P 1 

er b9 0 \J£ 


/. 


B„(X) 


calling it the "mirror point flow," or rate of change of mirror field per unit 
longitude. In absence of energy loss and sources, and taking into account (21) 
and (9), this quantity becomes 

v . = ? -7a 5 B^>-7 ?2 W (fcu> P 


in which 

<i> ~ p 2 =r 0 ^) < 27 > 

are mean changes of B and B 2 per unit longitude. 

Notice that there are three physically quite distinct contributions to the 
mirror point "flow.” The first term of the right hand of (26) represents a 
steady increase of B, (fi is positive), i.e., a steady lowering of mirror point alti- 
tude. This contribution comes from the first order, "streaming" term in B (21). 
In absence of dispersion (/3 2 = 0 ), is the only contribution to V m ; in that case, 

B = B (X) =Jv m dX is just the characteristic of the corresponding first order dif- 
ferential equation. 

The second and third terms on the right hand side of (26) arise in the dis- 
persion mechanisms. Their presence in (26) clearly shows that dispersion also 
contributes to a steady, average "flow" of mirror points along field lines. 

These two terms are governed by the gradients (along field lines) on the coeffi- 
cient/3 2 , and of U , respectively. The gradient of ft 2 is always positive, so that the 
second term in (26) always represents an upward motion of the particle's 
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mirror points (back-scattering from the denser atmosphere). If, on the other 
hand dU/^B < 0, the third term represents a lowering of the mirror point alti- 
tude. If, for instance, U/3 2 = const, along a line of force, the contribution from the 
two dispersion terms is zero: The number of mirror points thrown upwards by 
the gradient in/? 2 , exactly compensates the number of mirror points streaming 
downwards due to the gradient of U. The only mirror point flow which remains, 
in this case, is that of the first order term: V m = /?. 

Notice finally that 

C(B) =fui V m d 9 (28) 

represents the total number of electrons of our population crossing a constant-B 
ring per unit time. The integration is performed around the world along the 
B-L ring. It is easy to verify that (28) is independent of B, for a source-free 
diffusion. If V m is positive at all longitudes (lowering of mirror points), (28) 
represents the total number of electrons precipitating per unit time through a 
given constant-B ring. 

V. EXPRESSIONS FOR THE COEFFICIENTS 

In this section we shall obtain expressions for the various coefficients in- 
terventing in our fundamental equation (21). First of all, the half-bounce time 
T b is given in terms of the rectified path between mirror points (5): 

T b (B, ? ,E)=i (l + 2B^) (29) 

(I is usually computed in Earth Radii, so (29) must be multiplied byR e ). Ex- 
pression (29) has absolutely general validity, for any trapping field geometry. 


I- 22 



It offers the big advantage that it only involves computations of I. Results for 
the real geomagnetic field were cross-checked with numerical integration of the 
path S b and found to be coincident within high accuracy. For the angular drift 
velocity <p 0 , expressions (4) and (6) have to be evaluated numerically. 

The energy loss coefficient e is given by e = l/r b {e) where {e} is the aver- 
age energy loss per half-bounce: 


{e} 



BE 


ds 


B s 


/TT 


gjs) 

B 


(30) 


The origin of the field line arc length is taken at the equatorial point. ds/Vl -B'(syB 
is the element of trajectory of the electron at a point s where the field is B'(s). 

Using tables given by Berger and Seltzer, 21 we obtain the following very good 
approximation: 


where 



10 ' 2 ° N eff. ion. h ( E ) ( kev / cm ) 


> 


h(E) 


'2.61 + 349.0E* 0- 844 for 
4. 54 for 


E < 470 kev 
E > 470 kev 


J 


and 


N = 4. 17N» + N + 0.881 N . . + 

iN eff. ion. Argon oxygen nitrogen 


+ 0- 142N h,.,. S .o 


is the "effective" atmospheric number density for ionization loss of electrons. 
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As we are interested in low L shells only, we did not include the contribution 
from atmospheric ions and electrons. 

In order to evaluate {/3} , the average change of mirror point field B per 
half-bounce, we first have to perform the average of collisions over the isotropic 
azimuthal distribution of scattering angles. Following Welch, Kaufmann and 
Hess 11 and slightly changing their notation, we obtain 


2 = B(f, ~*) sin’S 

/Hs the change in mirror point field intensity B, when the electron's pitch angle 
scatters an amount <9 at a field position B' , averaged over all possible azimuthal 
angles of scattering. Now we have to find the average of sin 2 6 at the given field 
point B' . This will be given by 


where 



N i Vi 


r /2 , dcr. 

7j. = sin 2 e — j-d e 

j o dd 

is the average contribution from one atom of class i. N. are the number densi- 
ties of the different atmospheric constituents at the field point B' ; da./ dd 
are the differential cross sections for screened Coulomb scattering, for each 
constituent. Using Moli&re's expression for da/dO and following closely 
Welch, Kaufmann, and Hess, 11 we obtain 


L 


N; 


Vi = 4 X 10' 22 N 


ef f . 


s c a t t . 



+ 2 


— (cm" 1 ) 


( 33 ) 
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where 


N.„. ...... = S. 18 N a[gon + N OT>Een + 0.774 N nitIogen + 


+ 0.058 N. , . + 0.030 N k ... 

helium hydrogen 


(34) 


is the "effective" atmospheric number density for screened Coulomb scatter- 
ing of electrons. Again, the contribution of atmospheric ions and electrons 
is excluded. 

We finally have to perform the integration along the path of a particle from 
one mirror point to its conjugate: 


" B 

11 ^ 2 

Again following Welch, Kaufman, and Hess, 11 we obtain /3 : 

/S' 2 = 2B 2 (‘§> " sin2 6 


(35) 


Therefore: 


{/3 2 } = 2B 2 


L<b> ^ 'K? N ‘ 

" B 


(36) 


If we now call: 


S X (B, 9 ) = 


•(B) „ N eff.'scatt. ( s » 9) 


r KD> _B_ 
J-s(B) B '( S > 


A - b'(s>/b 


ds(cm' 2 ) 


S 2 (B,9) = 


■I 


S(B) N eff, scatt. O’ £ ) 

s(B) A - B'(s)/B 


ds (cm' 2 ) 
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S' 2 (B,<p) 


f S(B) N eff. ion (S.y) 

J-s(B) ^1 ~ B'(s)/B 


K(B, E) = 4 x 10‘ 22 B 



(C) C(E) = IQ" 20 h(E) (kev cm 2 ) 


we can summarize expressions (35), (36) and (30) in the form 

{/3} = K(S 1 - \ S 2 ) 

{/3 2 } - 2BK (Sj - S 2 ) ( 37 ) 

{€> =cs; 


In the region of very high atmospheric densities, as, say, below 200 km, 
the first two coefficients in (37) are related to each other. Suppose an effective 
density for scattering of the type 


N 


ef f. 


scatt. 


= N 0 (<P) exp 


( B - B a M) 

V AB(cp) / 


where B a (<p) i s the B value at the 100 Km level, and 


AB(<p) « B a (<p) 

With this expression, it can be shown that 
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and 


Sj^ 


S 2 - N 0 exp 


AB(9) 



Si 


AB(<p) 

S2 - Sl_ 2B- 


RB/3s is the average gradient of B along the field line, in the region of interest. 
Taking into account (55), the following relation between the coefficients (/ 3 ) and 
{/ 3 2 } holds: 

{/ 3 2 } = 2AB {/ 3 } (38) 

This is similar to a relation found by Walt and McDonald 12 for the coefficients 
of a time dependent (but longitude-independent) Fokker-Planck equation, set up 
in terms of mirror point altitudes, in the region of high atmospheric densities. 

Notice that with (38), the mirror point flow (26) at low altitudes becomes: 

V m - -AB^-| (4U) (39) 


Physically, the cancellation of the first term with the second term in (26) means 
that the upstream of mirror points due to diffusion from the denser atmosphere 
("backscattering" at lower altitudes) exactly compensates the downstream due to 
the first order term. What is left, is a mirror point streaming, entirely due to 
diffusion, and which may be upwards or downwards, according to the gradient of 
the actual electron distribution, along a line of force. A similar result was ob- 
tained numerically by MacDonald and Walt 13 . 

Notice at once that for a distribution function inversely proportional to the 
effective atmospheric density 
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which is close to a steady state solution of equation (21), expression (39) be- 
comes 

v.*J3 ( 40 ) 

at low altitudes. This is enough to predict qualitatively the behavior of the 
average "path" of the mirror point of an electron which drifts through the South 
American Anomaly: while the particle is still well West of the Anomaly, the 
mirror point will follow a constant-B ring. As the particle approaches the 
Anomaly , this ring will dip into higher atmospheric densities, and the mirror 
point will start "flowing" downwards (towards higher B values, because V m in 
(40) is always positive). If the particle was initially high enough in altitude (at 
low B), it will emerge again from the Anomaly, though with a higher mirror 
point B-value, eventually populating the B-L region which was voided in the 
Anomaly. After a sufficient number of drifts around the earth, the mirror point 
trajectory may start low enough on the West side, in order to succeed precipitating 
below the 100 km level in the Anomaly. Of course, during this whole process, 
the energy of a given electron decreases steadily. Electrons should mainly 
precipitate on the West side of the Anomaly. This purely qualitative analysis 
also shows that one should expect a diurnal variation of the electron distribution 
in the Anomaly, related to the diurnal variation of the atmosphere in that region. 

Notice finally, that according to this picture, replenishment should occur 
right at, or shortly "after," i.e., East of the Anomaly. Diffusion into the 
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"shadow" region far East of the Anomaly is negligible, for pure Coulomb scat- 
tering. We may conclude our discussion by pointing out that in the picture 
described above , there is a steady flow of particles out of the lower B region 
of a shell, towards higher B values, and from there, into the atmosphere at the 
Anomaly. In an equilibrium state, this loss must be exactly balanced by a con- 
tinuous injection from a source mainly effective in the lower B region. In ab- 
sence of such a source, i.e., in a non- equilibrium state, this flow and subse- 
quent loss should determine the lifetime of the trapped radiation in question. 

ACKNOWLEDGMENTS 

We are indebted to Drs. W. N. Hess, D. Stern, J. Herring and M. Liwshitz 
for valuable discussions and critical remarks. 

REFERENCES 

1. Yoshida S., G. H. Ludwig and J. A. Van Allen, "Distribution of Trapped 
Radiation in the Geomagnetic Field," Journal of Geophys. Res. 65, 807-813 
(1960) 

2. Vernov, S. N., E. V. Gorchakov, Yu. I. Logachev, V. E. Nesterov, N. F. 
Pisarenko, I. A. Savenko, A. E. Chudakov and P. I. Shavrin, "Investigations 
of Radiation During Flights of Satellites, Space Vehicles and Rockets, Suppl. 
J. Phys. Soc. of Japan 17, 162 (1962) 

3. Imhof, W. L., and R. V. Smith, "Longitudinal Variations of High Energy 
Electrons at Low Altitudes," J. Geophys. Res. 70, 569-577 (1965) 

4. Paulikas, G. A. and S. C. Freden, "Precipitation of Energetic Electrons 
into the Atmosphere," Journal of Geophys. Res. 69; 1239-1249 (1964) 

5. Mihalov, J. D., F. S. Mozer and R. S. White, "Artifically Injected Electrons 
at Low Altitude" Journal of Geophys. Res. 69, 4003-4013 (1964) 


1-29 


6. Freden, S. C. and G. A. Paulikas, "Trapped Protons at low Altitudes in the 
South Atlantic Magnetic Anomaly," Journal of Geophys. Res. 69, 1259-1269 
(1964) 

7. Heckman, H. H. and G. H. Nakano, "Direct Observations of Mirroring Pro- 
tons in the South Atlantic Anomaly," Space Research V, 329-342 (1965) 

8. Ghielmetti, H. S., N. Becerra, A. M. Godel, H. Heredia and J. G. Roederer, 
"Enhancement of the X-Ray Intensity at Balloon Altitudes in the South 
American Anomaly," Phys. Rev. Letters 12, 388-390 (1964) 

9. Williams, D. J. and J. W. Kohl, "Loss and Replenishment of Electrons at 
Middle Latitudes and High B-Values," Journal of Geophys. Res. 70, 4135- 
4150 (1965) 

10. MacDonald, W. M. and M. Walt, "Distribution Function of Magnetically 
Confined Electrons in a Scattering Atmosphere," Annals of Physics 15, 

44-62 (1961) 

11. Welch, J. A., R. L. Kaufmann and W. N. Hess, "Trapped Electron Time His- 
tories for L = 1.18 to L = 1.30," Journal of Geophys. Res. 68, 685-699 (1963) 

12. Walt, M. and W. M. MacDonald, "Diffusion of Electrons in the Van Allen Ra- 
diation Belt, I: Treatment of Particles with Mirroring Points at High Altitude," 
Journal of Geophys. Res. 67, 5013-5024 (1962) 

13. MacDonald, W. M. and M. Walt, "Diffusion of Electrons in the Van Allen 
Radiation Belt, II: Particles with Mirroring Points at Low Altitude," Journal 
of Geophys. Res. 67, 5025-5034 (1962) 

14. Anderson, A. D., G. E. Crane, W. E. Francis, L. L. Newkirk and M. Walt, 
"Theoretical Investigation of Geomagnetically Trapped Electrons from High 
Altitude Nuclear Explosions," Lockheed Missiles and Space Company Report 
LMSC 895355 (1964) 


1-30 


15. Mcllwaitt, C. E., "Coordinate for Mapping the Distributions of Magnetically 
Trapped Particles," Journal of Geophys. Res. 66. 3681-3691 (1961) 

16. Northrop, T. G. and E. Teller, "Stability of the Adiabatic Motion of Charged 
Particles in the Earth’s Field," Phys. Rev. 117, 215-225 (1960) 

17. Stone, E. C., "The Physical Significance and Application of L, B 0 and R 0 
to Geomagnetically Trapped Particles," Journal of Geophys. Res. 68, 
4157-4166 (1963) 

18. Roederer, J. G., "On the Adiabatic Motion of Energetic Particles in a Model 
Magnetosphere , " 

19. Northrop, T. G., "The Adiabatic Motion of Charged Particles," Interscience 
Publishers (1963) 

20. Hassitt, A., "The Drift Velocity of Trapped Particles," Journal of Geophys. 
Res. 70, 535-540 (1965) 

21. Berger, M. J. and S. M. Seltzer, NASA Report SP-3012 (1964) 


1-31 


LONGITUDE DEPENDENCE OF GEOMAGNETICALLY 


TRAPPED ELECTRONS 
PART II 


by 


Juan G. Roederer* 
Goddard Space Flight Center 
Greenbelt, Maryland 


Jasper A. Welch** 

Air Force Systems Command 
Los Angeles, California 

and 


James V. Herod*** 

Air Force Weapons Laboratory 
Kirtland AFB, New Mexico 


*National Academy of Sciences - National Research Council, Senior Post Doctoral Research 
Associate. Present address: Centro Nacional de Radiacion Cdsmica, Peru 272, Buenos 
Aires, Argentina 

**Lt. Colonel, United States Air Force 

***Captain, United States Army 


ABSTRACT 


The Fokker- Planck equation, discussed in Part I of this paper, was inte- 
grated numerically for a stationary case, for low L shells. A computer code 
was developed, which gives detailed information on particle shell geometry, and 
on longitude dependence of drift velocities and bounce periods. The coefficients 
representing longitudinal drift, ionization loss and multiple Coulomb scattering 
in the Fokker-Planck equation were computed numerically. 

The "eigenmode" electron distribution function for the longitude dependent 
problem was determined, and electron fluxes, energy spectra, B-dependence and 
mirror point trajectories were calculated for the area of the anomaly. The re- 
sults are in general agreement with satellite observations obtained in the vicinity 
of the anomaly. They do not, however, explain the observed replenishment which 
occurs at all other longitudes east of the anomaly. This latter effect must be 
caused by a pitch angle scattering mechanism other than multiple Coulomb in- 
teractions with atmospheric atoms. 
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LONGITUDE DEPENDENCE OF GEOMAGNETICALLY 
TRAPPED ELECTRONS 
PART II 


I. INTRODUCTION 

Pitch angle scattering is the main physical mechanism responsible for the 
shape of the equilibrium distribution of electrons in the inner radiation belt. 
These particles, if continuously injected into a given L-shell by a source, or 
transferred to it by radial diffusion, will gradually change their mirror points 
by such a pitch angle scattering mechanism, until a stationary state is achieved, 
in which the mirror point distribution represents an eigenstate of the corre- 
sponding mathematical diffusion problem. In this stationary state, injection 
rate equals the rate of removal, the latter being determined by the rate of 
transference of electrons to neighboring shells, and by the flow of mirror points 
into the dense part of the atmosphere. It is believed that for low L shells, 
radial transfer is much less efficient than atmospheric loss. 

In the inner magnetosphere there are at least two types of pitch angle 
scattering processes: multiple Coulomb scattering with atmospheric constitu- 
ents and resonant interactions with VLF electromagnetic or ion waves 1 ’ 2,3 ’ 4 . 
The first one must obviously be dominant at very low L values; the second one 
becomes more and more important as one goes to higher L shells. There are, 
however, reasons to expect that "electromagnetic" scattering is still effective 
at low L-shells, coexisting with Coulomb scattering. If this is the case, it 
should be difficult to determine experimentally how much of the scattering is of 
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non-Coulombic origin. The study of the longitude dependence of electron fluxes 
offers a possibility: Coulomb interactions are confined to a rather narrow band 
of longitudes in the South American Anomaly, and so must be the corresponding 
effect on the electron distribution. Electromagnetic interactions, instead, are 
not expected to have such a strong longitude dependence; the associated pitch 
angle diffusion should occur more or less uniformly over all longitudes. One 
possible way to obtain information about this electromagnetic scattering mech- 
anism at low L-shells is, therefore, to predict theoretically the Coulomb con- 
tribution to the longitude effect of the electron flux, and then to compare it with 
experimental results, blaming all observed differences on the electromagnetic 
scattering process. 

In Part I of this paper, the theoretical background for the study of the longi- 
tude dependence of a stationary electron distribution at low L-shells was given, 
under the assumption that multiple Coulomb scattering is the only interaction 
process present. In Part II we shall present results of numerical calculations 
for some particular L values, and compare them with the scarce experimental 
material so far available on this subject. 

H. FIELD GEOMETRY 

In order to solve the theoretical problem of longitude dependence of trapped 
electron distributions, it is first necessary to evaluate the coefficients inter- 
vening in the diffusion equation (21) of Part I. These coefficients, given in 
equations (37), contain integrals (S x , S 2 and S 2 ) of functions which vary many 
orders of magnitude over the domain of integration along the field line. There- 
fore, the first requirement for these calculations is an accurate description of 
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the geometry of magnetic shells, as well as the determination of the longitude 
dependence of all intervening kinetic variables, such as drift velocities and 
bounce periods. 

For particles mirroring at low altitudes in the South American Anomaly, 
changes of a few kilometers in altitude of magnetic shells may induce quite 
appreciable changes in the numerical results of theoretical diffusion calcula- 
tions. Considerable accuracy is therefore required. Of course, there will 
always be a natural limitation imposed by the degree of accuracy in the geo- 
magnetic field description. However, there are other more or less controllable 
sources of errors, both systematic and random, which are inherent to the 
numerical methods used in shell geometry and particle kinematics computations, 
and which must be reduced to a minimum. 

One of such sources is related to the question of shell splitting 5 . This shell 
splitting amounts to only a few kilometers at the equator, for L 1.2. Never- 
theless, such difference of a few kilometers may be of importance to mecha- 
nisms of atmospheric interactions, where it can compete with the scale height 
of atmospheric density. Special care must then be taken, when a particle shell 
is to be defined. For instance, if one is interested in the evolution of a group of 
particles mirroring far away from the equator, i.e. at high B-values, the field 
lines on which these particles will mirror must be defined by the proper B m , L 
(or B m , I) ring, and not by the equatorial ring of the corresponding L-shell. The 
differences which result from picking shells by their equatorial points instead 
of some higher B-value points, may be as much as 10-20 Km in the region of 
the Anomaly. The corresponding differences in longitude dependence of drift 
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velocities and bounce periods are, percentage wise, even higher (6-10% for 
L = 1.25); finally, the equatorial B-value of a shell picked at high B becomes 
appreciably longitude dependent (up to 6-7% for L = 1.25), and so does the 
equatorial pitch angle of a given particle mirroring at higher latitudes. 

Most of the longitudinal variations occur in the South American Anomaly, 
where, on the other hand, most of the atmospheric interactions take place. 
Unfortunately, this is the region where higher multipoles of the field expansion 
are felt most intensely, being therefore more exposed to errors of the particular 
field description used. 

A computer program was set up which meets all requirements of accuracy 
and speed 6 , giving detailed information on field line geometry for a given L-shell. 
In addition, it calculates the geometric factors intervening in drift velocities and 
bounce paths, using expressions (4) and (5) of Part I, respectively. As an example, 
the longitude dependence of the angular drift velocity (in arbitrary units) and that 
of the bounce path, for the L = 1.25 shell and two mirror point B-values, are 
shown in Figs. 1 and 2, respectively. When these curves are to be compared 
with previously calculated ones 7 ’ 8 , one has to bear in mind that longitude de- 
pendence of drift velocities and bounce paths very critically depend on how the 
particle shell was originally defined. The correct calculation requires picking 
the L-shell at precisely the B m points for which these quantities are to be 
determined. This is accomplished in the present code. 

III. CALCULATION OF THE COEFFICIENTS 

Once the field lines are defined on a given L-shell, the integrals for the co- 
efficients (37) have to be computed on each field line, and for different mirror 
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Fig. 1. Angular drift velocity as a function of the longitude of the equa- 
torial point, for particles mirroring at B = 0.18 and 0.22 gauss, for L = 1 .25. 



Fig. 2. Bounce path between conjugate mirror points, of electrons mirror- 
ing at B = 0.18 and 0.22 gauss, for L = 1.25. To obtain the half-bounce 
period, divide by the particle velocity in earth radii/sec. 
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point field B-values. These B-mesh points were defined at geodetic altitudes 
spaced roughly proportional to 1/2-1/5 of the local atmospheric scale height, 
at the low altitude end of each field line. The mesh was then conveniently com- 
pleted up to the equatorial point. On the other side of the equatorial point, the 
mesh was defined as the magnetically conjugate image of the original one. The 
B-mesh thus obtained had variable, longitude dependent cell size. In total, 100 
pairs of B-mesh points were defined on each field line. These field lines were 
spaced at longitude intervals of 5 or 3 degrees, decreasing to steps of 1° in the 
region of the anomaly. 

Atmospheric number densities intervemng in (32) and (34) depend on local 
time at altitudes above 150 Km. In order to take this into account, and to find 
its effect on the electron distribution, all computations were performed for two 
different local times in the region of the anomaly (0400 and 1400 LT). Each field 
line was assigned a mean local time according to its longitudinal position on the 
shell. The Harris-Priester atmosphere 9 was used; and two extreme cases of 
low and high solar activity were considered. For a given L- shell, the energy- 
independent parts of the coefficients were written on magnetic tape, for each of 
the four combinations night/day in the anomaly and low/high solar activity. A 
description of the computer code will be published elsewhere. 

IV. NUMERICAL SOLUTION OF THE DIFFUSION EQUATION 

As a first step, mainly to check the coefficients calculated in the previous 
section, the time dependent equation (19) was solved using coefficients averaged 
over local time and longitude as prescribed in expression (20). Results are in 
general agreement with previous calculations 10 , in particular, for the decay of 
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the Starfish injection 11 . The disagreement with experimental results for higher 
B-values, found by Welch et al. 10 and by Walt and Newkirk 11 , still subsists in 
our calculations (theoretical calculations always give much faster decays at high 
B-values, than experimentally observed). 

We now turn to the numerical solution of equation (21), for a stationary state. 
In order to integrate equation (21), the source term Q must be known. This, un- 
fortunately, is not possible for the time being. However, injection during inter- 
vals of time comparable to a typical drift period must be expected to be rather 
small at low L-shells. We shall therefore neglect the contribution from Q, 
always keeping in mind that with this assumption, a true steady state can never 
be achieved in the numerical calculations. The error committed should be small 
in comparison with other inaccuracies (atmospheric density, magnetic field, etc.). 

To solve numerically eq. (21) under these conditions, we have to start with 
a given initial distribution U at an initial longitude <p 4 . Proceeding once around 
the world with the numerical integration, until the starting longitude is again 
reached, one ends up with a distribution which, in general, will be different from 
the initial one. In particular, in absence of a source term Q, the integral flux 
of particles will have diminished due to the action of the sink at low altitudes. 

We shall arbitrarily call an "eigenmode" of eq. (21), a distribution for which 
this "one-time-around" difference is minimum (for a true eigenmode, this dif- 
ference is zero). 

The point is, then, to start at <p. with a distribution which already is as 
close as possible to an eigenmode and to proceed integrating around the earth 
a sufficient number of times until the "one-time-around" difference have be- 
come small enough. 
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For practical purposes it is not necessary to integrate every time around 
the whole world. It will be sufficient to consider a longitude interval, centered 
at the anomaly, with a width determined by the range of B-values of interest. 

For instance, for the L = 1.25 shell, it is sufficient to consider the longitude 
interval between approximately -70° and -18°, if one wants to explore the region 
B 1 0.270 gauss. Outside this longitude interval, the effect of multiple Coulomb 
scattering is negligible, due to the high altitude of the B-L rings under 
consideration. 

The procedure to obtain the stationary, longitude dependent electron dis- 
tribution is sketched in Fig. 3, for the case L = 1.25. We have used the solution 
of equation (19) obtained for long periods of time after injection, as the initial 
distribution for the longitude dependent equation (21), injected right at the center 
of the anomaly (cp.). The output distribution obtained at cp = -20° is reinjected at -70°, 
and this process is then repeated several times, until the percentage change of U 
per integration cycle (AU= U_ 20 - U_ 70 ) reaches a constant value. The result- 
ing distribution U is then adopted as "eigenmode" of the longitude dependent 
problem. In the real case, the small differences AU (<0) must be cancelled by 
the longitude integral of the injection rate jQdcp. Numerical results have shown 
that, in general, the distribution U which builds up in the region of replenish- 
ment east of the Anomaly, is only very weakly dependent on the initial distribu- 
tion injected at cp. . In other words, the eigenmode is achieved very fast in this 
high-B region. 

The computer program for the integration of eq. (21) was designed in such 
a way as to perform the reinjections as many times as wanted. In this program, 
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the initial distribution is read in, and the energy dependent parts of the coeffi- 
cients are computed. The field-geometric part of these coefficients is then read in 
from magnetic tape, for the first pair of field lines, and the associated difference 
equations are solved. This procedure is continued until a prefixed longitude east 
of the anomaly is reached, at which integration stops. The result is reinjected 
at another fixed longitude west of the anomaly, and the whole process is repeated. 

In order to better interpret the physical implications of the procedure to in- 
tegrate eq. (21), outlined above, we go back to the general equation (16). If we 
suppose a sudden injection at a time t i and a longitude cp if of the form 

Q = U, S(t-t i )S(cp-qp i ) 

and integrate equation (16) from t. - 8t to infinity, we have, taking into account 
(17): 


T b^0 


£ {■%) = - 4 > - k > <41) 


where 


■1 


W(B, cp, E) = U(B,<p,E, t)dt 


(42) 


According to (3) and (42), the quantity 

B 2 




represents the total number of electrons of energy between E and E + dE which 
have passed once in their longitudinal drift through the longitude cp , per unit area 
perpendicular to their pitch angle. This clearly defines the physical meaning of W. 
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Notice now that eq. (41) is formally identical to the longitude dependent 
equation (21), in absence of a source term Q. In conclusion, the solution of 
eq. (21) for Q = 0, and for a given initial distribution U injected at <P A , physi- 
cally represents, at each longitude, the total number of electrons which would 
have passed through that longitude on their first trip around the earth after 
injection (it does not contain any information about electrons going through on 
their second, third, etc., time around, nor does it say what the distribution was 
at any given time t). 

V. RESULTS 

Before analyzing general results, we turn to the discussion of the mirror 
point flow, given in section IV of Part I. Once the solution of eq. (21) has been 
obtained, one can integrate (25) in order to find the "mirror point trajectories" 

B m = B m (<P). Fig. 4 shows the result for 600 kev electrons, for L = 1.25 (0400 LT 
in the anomaly, low solar activity). All qualitative predictions, given in the last 
section of Part I are clearly confirmed. Notice the "window" in B which feeds 
the B-region originally wiped out of particles in the anomaly, and watch how 
higher mirror points are gradually lowered in altitude, as they drift through the 
anomaly. Constant altitude curves are reproduced in order to show at which 
altitudes most of the mirror point flow occurs. It can be seen that the slope (V m ) 
of the mirror point trajectories is nearly constant along constant altitude curves. 

Figs. 5 and 6 show general results of the integration of eq. (21). The com- 
puted, directional flux j Q of 300 kev electrons is plotted for different L-values 
as a function of longitude, normalized for different B-values outside the anomaly. 
Notice the following features: (a) Replenishment via multiple Coulomb scattering 
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Fig. 4. Mirror point trajectories of 600 kev electrons (average energy of electrons in the region 
of replenishment). Constant altitude curves are shown. Notice the narrow ‘‘window of mirror 
point field intensities around B = 0.218 gauss, west of the anomaly, which feeds t e region 
of replenishment on the east side. 
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occurs within 20-25° east of the point of emergence from the 100 Km level of 
the corresponding B-L ring, (b) The process is not symmetric with respect to 
the anomaly: scattering out (loss) of a B-interval west of the anomaly is faster, 
than scattering in (replenishment) on the east side. This asymmetry is more 
pronounced than the natural asymmetry of B-L rings, (c) At lower L values, 
for which B-L rings submerge into, and emerge from, the dense atmosphere in 
much steeper fashion, a "piling up" of particles occurs just prior to the decrease, 
and just after replenishment. This piling up shows up, but less pronounced, in 
the omnidirectional fluxes. 

Next, we reproduce in Figs. 7 and 8 the energy spectra of electrons mirror- 
ing at different, high B-values, as they appear in the eigenmode distribution, 
outside the anomaly. An inspection of these figures reveals that replenishment 
is more efficient at lower energies. For electrons of £ 2 Mev, replenishment 
by multiple Coulomb scattering east of the anomaly is negligible. Whatever high 
energy electron flux is found east of the anomaly in the high B region, must have 
diffused in by means of a process other than multiple Coulomb scattering. 

Dependence with B of the omnidirectional flux is very steep; some typical 
curves are shown in Fig. 9. Finally, the omnidirectional flux of >300 kev elec- 
trons at a fixed altitude of 300 km is plotted as a function of longitude for two L 
values in Fig. 10. 

In analyzing these results, one has to keep in mind that in the region of 
replenishment (high B), the distribution is almost completely independent of the 
originally injected spectrum, and of all what happens near the equator, if multiple 
Coulomb scattering is the only process involved. 
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Fig. 7. Differential energy spectra of electrons mirroring at high B-values west 

of the anomaly, for L = 1.18. 
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Fig. 10 Variation of the omnidirectional flux of >300 kev electrons along a 
constant altitude level (300 km), for L = 1.18 and 1.25. 
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VI. COMPARISON WITH EXPERIMENT 

Detailed experimental data on longitude dependence of electron fluxes at 
low L shells are not very abundant. Only two papers deal explicitly with this 
subject; both describe results obtained with polar orbiting satellites. In the 
first one, Imhof and Smith 12 present results on B and E dependence of omni- 
directional electron fluxes, grouped into two longitude intervals, called "east" 
(-30° < op < + 10°) and "west" (-90° < cp < - 60°) of the anomaly, for different 
L-shells (L < 1.7). In the second paper, Williams and Kohl 13 study electron 
fluxes for L > 1.6 at three particular longitudes: cp = 0°, 211° and 282°. 

At the present state of experimental information, only more or less quali- 
tative comparisons with our calculations are possible. Let us consider the 
L = 1.25 shell, and the measured B-dependence of the omnidirectional flux 
(second graph in Fig. 4 of ref. 12). According to our Fig. 6, the electron dis- 
tribution in the longitude interval -30° < cp < + 10° ("east" of the anomaly) 
should be almost completely "Coulomb-replenished" for the L = 1.25 shell. 

The computed curve for > 300 kev electrons in Fig. 9 la in very good agreement 
with the experimental B-dependence, except for high B-values (£0.30 gauss, 
which for field geometry reasons can have been sampled by the satellite only 
near <p = 10°). 

Imhof and Smith's results for -90° < qp < - 60° ("west" of the anomaly) 
clearly indicate the action of an additional pitch angle scattering mechanism. 
The authors suggest large-angle Coulomb scattering, by means of which equa- 
torial particles would be scattered into high-B mirror point trajectories. We 
wish to point out that such a process, also atmosphere-controlled, would lead 
for low L shells to a replenishment process almost as longitude dependent as 
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multiple scattering, being most effective in the region of the anomaly. Experi- 
mental results on electron fluxes and spectra in the region cp ~ 120° are urgently 
needed for more conclusive studies. 

We now turn to the results obtained by Williams and Kohl. Unfortunately, 
they refer to observations made at higher L-values than those under considera- 
tion in our theoretical and numerical analysis. One of the general conclusions 
of these authors, regarding >280 kev electrons, is that electrons mirroring at 
B-values which remain above 350 km in altitude over the South Atlantic anomaly 
are not disturbed by the anomaly; about 90% of those that reach the 100- to 250 
Km region over the anomaly are lost, and all those electrons reaching an altitude 
of <100 km over the anomaly are lost. An inspection of the mirror point tra- 
jectories reproduced in Fig. 4 broadly confirm this. 

Williams and Kohl present results picturing the "transmission" of electrons 
through the anomaly for different L shells. They define the "transmission" as 
the ratio of the counting rates obtained at cp = 0° and cp = 282°, respectively, for 
a given pair of B,L values. Considering the L range explored, the cp = 0° count- 
ing rates really correspond to the "inside" rather than to the east side of the 
anomaly, (see their Fig. 5). This means that east of their cp = 0° point, Coulomb 
replenishment should still be going on for some 20° to 30° (remember that for 
higher L values, B-L rings emerge much slower out of low altitudes). True 
transmission factors may then be considerably higher than those quoted by 
Williams and Kohl. This is why a comparison of their results with our curves 
in Figs. 5 and 6 is not too meaningful. 
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Comparison of their <p = 282° and cp = 211° data again clearly indicates the 
action of an additional, non-Coulombic pitch angle scattering mechanism, acting 
well east of the anomaly. The authors suggest resonant interactions with electro- 
magnetic VLF waves. A most puzzling result is that in the high-B region, they 
find particle fluxes increasing by a factor 2-3 in the longitude interval between 
211° and 282° (i.e. only 71° wide), for L ~ 2.0. This result could never be ex- 
plained by a pitch angle scattering mechanism acting with uniform efficiency 
over all longitudes. For the time being, it is difficult to understand why the 
additional scattering mechanism should be so much more effective in the 211- 
282° longitude interval. 

Electron spectra obtained in our calculations (Figs. 7,8) are very soft in the 
region of replenishment. This is in qualitative agreement with observations by 
Paulikas and Freden 14 . 

The curves given in Fig. 10 are in quite poor agreement with similar plots 
extracted from isoflux contours shown by Mihalov et al. 15 and by Vernov et al. 16 . 
The experimentally observed longitude dependence of the omnidirectional flux at 
a fixed altitude for a given L-shell, is much broader than the narrow "peaks” 
shown in Fig. 10. However, one has to keep in mind that the Mihalov et al. data 
were obtained for artifically injected electrons (much harder spectrum), and 
that Vernov et al. present GM counting rates, i.e. a superposition of electron 
and proton fluxes. Notice that wider peaks in a representation like Fig. 10 are 
equivalent to a less steep B-dependence of the flux along a field line (Fig. 9). 
This, again, points to the action of a non-Coulombic scattering mechanism at all 
longitudes. 
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VII. CONCLUSIONS 


Coulomb scattering is certainly the most important pitch angle diffusion 
mechanism for low L- shells, but not the only one. Presently available experi- 
mental material is not detailed enough in order to reach more quantitative con- 
clusions about the extra scattering mechanism. 

It would be very useful to have more information on omni- and unidirectional 
electron fluxes, their energy spectra and B-dcpendence, at the following longi- 
tudinal positions for L < 1.3: at the "center" of the anomaly, at 15-20° east and 
west, respectively, and at 60-90° longitude intervals, in the rest. This informa- 
tion should be grouped into two broad UT intervals (corresponding to day and 
night in the anomaly). When this information is available, it would be much 
simpler to "filter out" the theoretically predicted multiple Coulomb scattering 
contribution, and to study the behaviour of the remaining effect. 
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FIGURE CAPTIONS 


Figure 1 . Angular drift velocity as a function of the longitude of the equatorial 
point, for particles mirroring at B = 0.18 and 0.22 gauss, for L = 1.25. 

Figure 2. Bounce path between conjugate mirror points, of electrons mirroring 
at B = 0.18 and 0.22 gauss, for L = 1.25. To obtain the half-bounce period, 
divide by the particle velocity in earth radii/sec. 

Figure 3. Sketch to illustrate the procedure of numerical integration of 
equation (21) (Part I). 

Figure 4. Mirror point trajectories of 600 kev electrons (average energy of 
electrons in the region of replenishment). Constant altitude curves are 
shown. Notice the narrow "window" of mirror point field intensities around 
B = 0.218 gauss, west of the anomaly, which "feeds" the region of replen- 
ishment on the east side. 

Figure 5. Longitude dependence of the directional flux j of 300 kev electrons 
mirroring at different B-values, for L = 1.18. 

Figure 6. Longitude dependence of the directional flux j of 300 kev electrons 
mirroring at different B-values, for L = 1.25. 

Figure 7. Differential energy spectra of electrons mirroring at high B-values 
west of the anomaly, for L = 1.18. 





Figure 8. Differential energy spectra of electrons mirroring at high B-values 
west of the anomaly, for L = 1.25. 

Figure 9. B-dependence of the omnidirectional flux of electrons, for various 
energies. 

Figure 10. Variation of the omnidirectional flux of > 300 kev electrons along a 
constant altitude level (300 km), for L = 1.18 and 1.25. 
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